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Abstract 

We discuss a possible new manifestation of the formation of an initial 
nonequilibrium flux tube stage in high energy heavy ion collisions. It is 
shown that a strong asymmetry in particle distributions around axes of 
hard-transverse jets takes place due to a large difference between longitu- 
dinal and transverse forces acting on the hardly scattered quarks or gluons 
crossing the flux tubes. 

Jet formation and evolution in hot-dense hadron matter has been 
a subject of recent papers. This problem is considered in many 
studies: interaction of the leading partons of di-jets with the ex- 
cited matter produced in ultrarelativistic nucleus-nucleus collisions 
leads to energy-momentum transfer to the hot and dense medium, 
which may manifest itself as jet quenching [|I], @], a broadening of 
the acoplanarity distribution || f|, |5| and possible appearance of 
dynamic instabilities in quark-gluon plasma (QGP) |6|, |7|, ^ which 
leads to characteristic pion and photon emission. Also the quantum- 
chromodynamics (QCD) flux tube dynamics at finite temperature 
can lead to broader pr distribution and increased strangeness yield 

In this paper we suggest a new signature of the formation of an 



initial flux tube like stage |T0[ . It is shown that presence of the flux 
tube stage in collisions leads to a modification of distributions of 
secondary particles in hard transverse jets. In contrast to a single 
jet case, in which the distributions are axially symmetrical relative 
to the jet axis, the modification appears as the symmetry violation 
(Fig. p. This phenomenon should take place due to the large dif- 
ference between the longitudinal and transverse forces acting on the 
hardly scattered and almost transversely moved quarks or gluons 
which crossed the flux tubes and formed the jets. So QCD hard jets 
can serve as probes of the hot and dense matter in the flux tube 
stage. 

Let us assume that formation of the multi flux tube state in 
nucleus-nucleus collisions takes place at the initial stage: during 
the collision of two heavy nuclei a system of flux tubes begin to 
form in the volume of the overlap of the colliding nuclei. At a later 
stage, if the conditions reached in relativistic heavy ion collisions 
are adequate, colorless QGP may be formed via fusion of the tubes 
and the color would no longer be confined to hadronic or flux tube 
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Figure 1 : An illustration for the asymmetry of the secondary particle distribu- 
tion in the hard jets and the initial nonequilibrium flux tube stage. The cone 
of dashed lines corresponds to the jet in the case of no flux tube stage and the 
cone of full lines corresponds to the asymmetric jet. Black circles denote quarks 
(triplet color representation) and white circles denote antiquarks or diquarks 
(antitriplet representation) . 



3 



dimensions. 

This flux tube initial state, if created, immediately begins its 
complicated evolution to ordinary hadronic matter — directly or 
passing through the QGP state. The question is, which probes can 
provide direct information about this flux tube crossover form of 
the matter? Considering the longitudinal flux-tube-like field fluc- 
tuations oriented along the reaction axis which cannot be realized 
in a pure hadronic case or in a case of equilibrium QGP, we suggest 
a new backgroundless signature of the formation of the initial stage 
in contrast to most known experimental phenomena which do not 
give clear information about the pre-equilibrium stage. 

Let's consider the flux tube stage for the central nucleus-nucleus 
collision (Fig. [1]). The tube-like fluctuations of the color fields are 
oriented along the collision axis. This anisotropy of the field fluctu- 
ations affects the secondary particle distributions. 

For a first estimation in a frame of the classic approach, consider 
a hardly scattered parton which is moving in a direction transverse 
to the reaction axis. This parton may create a flux tube which 
may decay to a hard jet. An interaction of this parton with QCD 
fields in the longitudinal flux tubes affects the parton's momentum 
distribution which determines the distribution of secondary particles 
in the jet relative to the jet axis. In particular, bremsstrahlung gluon 
radiation takes place in a tangential direction by the quark crossing 
the strong gluon field in each tube. This is the main dynamical 
reason for an asymmetry: the mean square projection of momentum 
of secondary particles in the jet along the reaction axis (Fig. [I], 
axis z) becomes larger than in the transverse direction (axis x). 
Because the fluctuations of the fields are isotropic in a vacuum, 
hadron matter, or QGP states, the distribution of the jet particles 
in these cases is axially symmetrical relative to the jet axis. 

Our estimations show that this asymmetry can be visible in 
nucleus-nucleus collisions if the flux tube stage takes place. The 
value of the asymmetry can be characterized by the ratio R a = 
<p^> 1 ^ 2 /<p^> 1//2 . Let's estimate the value of the asymmetry of the 
parton distribution relative to axis of its initial direction of move- 
ment (immediately after the hard scattering). 

R a depends on energy of colliding nuclei because of energy de- 
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pendence of the mean number of flux tubes and mean square 
transverse momentum of particles in a jet relative to jet axis. 
Consider e + e~-annihilation to estimate the energy dependence of 
<P 2 > 1 ^ 2 — <P 2 Z > 1 ■ In this case charge particle multiplicity can be 
parametrized as fill 



M ch (W) = a + b\nW 2 + c\n"W\ (1) 

where a = 3.33, b = —0.40, c = 0.26 and W is invariant mass of 
e + e~ pair. The total particle multiplicity can be roughly estimated 
as M(W) = 3 /2M c]l (W), and the value of <p 2 > 1//2 can be expressed 
from sphericity S ||TT| 

c _ 3 J2i(pj X + Pjz) / Q \ 

6 " 2 T,iP 2 ' (2) 

and the asymptotic dependence of W [[L2| 

S oc 0.8W^ 1/2 . (3) 

If the jet is not too wide (<p 2 > + <p 2 > <C <Py> — <P 2 >, see 
Fig. 0), then W « M<p 2 > x / 2 and 

W=M^V=^. (4) 

From (f|) we get 



<pl> 1/2 = \\<pI + p> 



w 3 / 4 

E(pi+pL)^o.52- 



\ 2M Y 

(5) 

For SpS, RHIC and LHC energies we obtain <#> 1/2 = 0.26, 0.44 
and 1.8 GeV/c respectively. 

Now let's estimate the additional momentum of the hard trans- 
verse quark crossed the longitudinal flux tubes. Consider the sim- 
plest case of static cylindrical flux tubes in ground states with a 
constant uniform chromoelectric field pointed along the axis of sym- 
metry (the reaction axis). 

Then consider a quark (qi) passing through the flux tube trans- 
verse to the tube axis. Static tube approximation can be used for 
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Figure 2: An illustration of the calculation of quark momentum and for discus- 
sion of the interaction between a hard scattered quark and flux tubes. 



the piece of the tube (crossed by the quark) which is almost at rest 
in the equal velocity frame. Because in static tubes chromomag- 
netic fields vanish, and soft interactions of the hard quark qi with 
its own (transverse) tube (Fig. ^]) can be neglected for a study of the 
transverse dynamics, the force acting on the quark can be written 
as 

F = |a?E b , (6) 

where g is the QCD coupling constant, A a are the Gell-Mann ma- 
trices, the lower index 1 means that the matrix acts on qi and E a 
is the chromoelectric field. For simplicity, the flux tube lengths are 
assumed to be much greater than their radii, and the flux tube radii 
to be small enough to neglect the overlap of the tubes. For the 
quark induced case the tube ends have triplet color representations. 
The quark qi moves with almost light velocity, therefore, in the thin 
tube case, the end quarks have no time to change their color states 
when qi crosses the flux tube. So the values of E a are distributed 
among the tubes in accordance with their color states and can be 
considered time independent. 

Find the change of the mean transverse momentum of the hard 
quark after crossing one flux tube with an arbitrary distribution of 



G 



the gluonic field over the tube section 



E a = E a (x,y). (7) 

Let axes x and y be the coordinate system in the section where qi 
crosses the tube, let axis y be pointed along qis direction of motion 
and let £ be the distance between axis y and line of q\ movement. 
Then q\ will cross the tube surface in points (£, yi(£)) an d (£, 2/2(0)- 
The quark momentum change may be written as 

2/2 (0 2/2 (£) 

Ap(0 = /F(f,t)dfc= | F(^ 2 /)d 2 /= | ^AJE'&j/Jdy. (8) 

yi(0 2/1 (0 

Here and below the light velocity (c) is assumed to be c = 1. Average 
Ap over all £ (ptube is the flux tube radius): 

Ptubc 1/2 (0 

AP = / « / <7-A?E"(£,y)^ = 

J «<S-Ptubo 2/1 (?) 

— Ptube 

1 J g l -XlE\x,y)dS. (9) 



2p. 



tube 



Here S 1 is a transverse flux tube section. Applying the Gauss theo- 
rem (see Fig. |) we get: 



a p = ^4 A i/ Earfs =^4 A i/ 

z Ptube z £ z Ptube z ^ 



div E a dV 



z -^A;^ = i-^-A;AS, (io) 



2ptubc 2 2 8p tu bc 

where g^A^ is the "color charge" of the quark q 2 and z is the unit 
vector in the longitudinal tube direction. This result does not de- 
pend on the field distribution over the tube section. 

For estimation of the asymmetry of distribution of particles in 
the jet mentioned above we need to know <Ap 2 >: 

<A P 2 > = ^<(^A«) 2 >, (11) 
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Here the index c means averaging over the color states of q\ and q 2 . 
To calculate <(A"A2) 2 > C let's consider matrices A" and X 2 as vectors 
in eight- dimensional space. Choose the axis x$ to be directed along 
the vector A" (its distribution is assumed to be isotropic). Take 
into account that the mean square of the scalar product (AfAf) 2 
averaged over all A| directions and the mean value of the square 
of the projection of Af onto the direction of A" (i.e. <(A 8 ) 2 > C ) 
multiplied by the absolute value of A" squared are equal to each 
other. This takes place because the color states of q\ and q 2 are 
assumed to be independent. Taking into account that all mean 
square components of the vector A a are equal to each other we get: 

<A 2 > C = <A a A a > c = ]T<(A a ) 2 > c = 8<(A 8 ) 2 > c , 

a=l 

<(A 8 ) 2 > C = ^<A 2 > c = iA 2 . (12) 

For A 2 = f : 

<m2>c 4(i) 2 - (13) 

Let Ap = \/ <Ap 2 > c . Then after crossing one flux tube 

2 

where a s = |- is the QCD fine structure constant. 

For a tube radius ~ 0.3 Fm [O we get for QCD running cou- 



pling constant a s( AQC ^ Pt b ) = 0.3 (in 3-loop approximation) and 
Ap w 0.6 GeV. In the cases of SpS and RHIC (see Table || be- 
low) this value is greater than or close to <p 2 x > 1 ^ 2 and asymmetry 
R a ~ 1.7-2.5 is larger than in the case of LHC (R a ~ 1.05). 

This asymmetry is expected to be visible for hadron-hadron col- 
lisions, and increases together with the atomic number A of the 
colliding nuclei because the quark crosses more and more tubes. In 
this paper we consider the jet asymmetry in the large A case for the 
reaction of 238 U+ 238 U. 

Now let's estimate the average number of flux tubes crossed 
by the quark in its traveling across the reaction zone in central 
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nucleus-nucleus collisions. For kinematical reasons in addition to 
the quark qi there is a quark q[ hardly collided with q\ and moving 
in the opposite direction from q\. The mean distance travelled by 
each quark in the nuclear reaction zone with uniform distribution of 
the flux tubes is then 

<'> = £ < 15 > 

Here R is the reaction zone radius (in the case of the central in- 
teractions of nuclei it is almost equal to the radius of the lightest 
nucleus) . 

We assume that the mean number of flux tubes is about 3/2 
for each nucleon because for each nucleon there are four possible 
numbers 0, 1, 2, 3 of the tubes encountered with almost equal prob- 
ability. This rough estimation is in good agreement with FRITIOF 
]14| calculations. 



Consider a central collision of two identical nuclei with atomic 
number A. In this case the number of flux tubes per unit area of 
the transverse section is 

n= ^=2^r (16) 

Here R = tqA 1 ^ 3 is the nucleus radius. Then the average number of 
flux tubes crossed by q± will be 

8 ^2/3 

N = 2p tnhe <l>n = — — p tube . (17) 

7T 7™ o 

For 238 U we get iV « 8. 

Let's find the average transverse momentum of quark qi after it 
crosses N flux tubes. After crossing m flux tubes 

<AP z 2 m > = <{AP zm _ 1 + Ap z ) 2 > = 

= <AP^ 1 > + <A Pz 2 > + 2<AP* m _ 1 A Pz >.(18) 

Here AP zm _ x is the transverse momentum of q\ after it crossed m — 1 
flux tubes. The last term vanishes because AP zm _i and Ap z may 
be pointed in any direction along the tube axis independently. So, 
we get 

<AP 2 m > = <AP Z ^ 1 > + <A Pz 2 >. (19) 
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And 

<AP> = J<AP 2 N > ~ Jn(A p ) 2 = VNAp (20) 



AD 8A 2 / 3 Ana s (p tuhc ) 2-kA 1 / 3 a s (p tube 

<^-r > — 4 / — p tube 



vr 2 r 3\/2p tube v 7 ^) V^tubc 

« 1.7GeV/c. (21) 

This estimation is accurate for the case of non-overlapping flux 
tubes. If the flux tube radii are not too large then this assumption 
is quite reasonable. 

Let's estimate the asymmetry of particle distribution in a jet 
characterized by R a = <p 2 > 1 / 2 / ' <p 2 x > x l 2 . The value of <p 2 > 1 / 2 
is determined by average transverse momentum of the hard quark 
<AP> and by typical transverse momentum of particles in the 
jet <pl> 1 ' 2 : 



Ra 



<AP> 2 



^1/2 ^<AP> 2 + <p 2 x > 

<pl> 1/2 ~ <p> 1/2 ~ Y ' <pI> 

(22) 



16 A 2 / 3 a 2 



\ 3 r p tubc <p2>' 

The values of asymmetry R a at various energies of colliding nuclei 
are shown in Table ^| (see below). One can see that in the cases of 
SpS and RHIC R a ^> 1 and so large asymmetry can be measured 
experimentally. In the case of LHC the asymmetry is not very large 
but may be observed in high statistics experiments. 

It should be noted that our calculation of quark interaction with 
flux tubes is not complete yet. For the first step we concentrated 
on the subject of phenomenon and did not take into account the 
quark's interaction with fields in its own flux tube or the interactions 
between tubes. We consider only quark jets here although we would 
expect the more noticeable effect for gluon-induced jets. The flux 
tubes are assumed to be non-overlapping. In Table [1] we present 
the results of our estimation (based on event generator FRITIOF) 
of the mean numbers of essentially overlapped tubes for different 
tube radii for the reaction of 238 [f + 238 U. We suppose here that 
the tubes are essentially overlapped when the distance between their 
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Table 1: Mean numbers of n— fold overlapped flux tubes for different tube radii 
(n = corresponds to the not overlapped single flux tube). 





Tube radius (Fm) 


n 


0.1 


0.2 


0.3 


0.4 





330.9 


261.6 


175.9 


102.7 


1 


25.4 


82.0 


123.9 


120.9 


2 


0.7 


11.6 


43.5 


81.3 


3 


0.0 


1.6 


10.5 


35.3 


4 


0.0 


0.2 


2.3 


12.0 



axes is less than the tube radius. The resulting value (^) takes into 
account only the most important contributions and needs in more 
precise calculations which we plan to make in our following works. 
The small value of tube radius 0.2-0.3 Fm are of great interest (in 
connection with percolation [T5| and lattice calculations One 
can see from the Table [l] that in this case the most of tubes are not 
overlapped. 

It should be noted that a jet's properties can be studied in the 
experiment if its angular size is sufficiently large to cover more than 
one detector module and if there is no dominance of essentially over- 
lapped jets. 

To estimate the angular size of the jets let's consider a jet in 
the equal velocity frame where the partons have been elastically 
scattered off each other at the angle of 90 degrees. Each parton 
carries about x — 1/3 of the total nucleon momentum p (sea quarks 
and gluons do not change strongly this value). The angular size of 
a jet is determined by (Fig. ||) 

. e <pI> 1/2 ^<pI> 1/2 m(w) 2 1/2 

where <p 2 x > 1 / 2 is shown in (|5]) and the factor 2 in the denominator 
appears because M is the multiplicity of the two jets propagating 
in opposite directions. In Table |2| we present values of po, <p 2 > x ' 2 
<p 2 > 1 / 2 , R a and for SpS, RHIC and LHC energies. One can see 
that in all cases the angular size of the jets is large enough to be 
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(xi = x 2 = 1/3) 




detected experimentally. 

We have yet no estimations of the numbers of overlapping jets at 
different transverse momenta. But in any case there are transverse 
jets with high momentum which are rare enough to be single. 

Other corrections may arise from the decay of the longitudinal 
tubes on fragments due to their rotation in evolution of the flux-tube 
stage. The largest asymmetry takes place in the case when the quark 
crosses the flux tubes transversely. Because the angular velocity of 
tube rotation is greater in case of short tube fragments (the mean 
square transverse momentum of fragment ends does not depend on 
the fragment lengths), events in which the long tube fragments dom- 
inate are more favorable for studying the asymmetry. In particular 
such events may be selected by using the fact that longer tubes ac- 
cumulate more kinetic energy in their tension and their fragments 
after decay should be concentrated in the middle rapidity range. 
We plan to present more detailed calculations of this effect in our 
future works. In any case the large value of asymmetry (~ 670%) 
means that this asymmetry can be observed in experiments even if 
the corrections would be taken into account. 

Another reason for decreasing of the asymmetry is the percolation 



effect. As shown in [y| even in central Pb-Pb collisions string 
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Table 2: Nucleon momentum pg, invariant mass of hardly scattered partons W, 
particle multiplicity of di-jet M, typical scale of the particle transverse momenta 
in jets in the beam axis direction (<p^> 1 / 2 ) and in direction perpendicular to 
the jet and beam axes (Kp 2 ^ 1 / 2 ), value of jet asymmetry R a and jet angular 
size 9 for different colliders. 



Collider 


SPS 


RHIC 


LHC 


Po, GeV/c 


10 


100 


3150 


W, GeV 


20/3 


200/3 


6300/3 


M 


8.33 


27.5 


87.1 


<p 2 x > 1/2 , GeV/c 


0.26 


0.44 


1.8 


<P 2 Z > 1/2 , Gev/c 


1.7 


1.7 


2.5 




6.7 


4.0 


1.4 


o,° 


38 


21 


8.7 



density is so high that practically all tubes are fused. It means 
that there is no jet asymmetry. But it is not yet known if this 
picture is real or not. Observing any asymmetry or its absence in 
experiments may give information about the way of evolution of 
heavy ion collisions. 

In conclusion it should be noted that a new possible manifestation 
of the formation of the initial nonequilibrium flux tube stage for the 
QGP transition in high energy heavy ion collisions is suggested: 
the longitudinal anisotropy in the flux tube fluctuations of color 
electric fields formed in the pre-equilibrium (flux tube) stage leads 
to a strong azimuthal asymmetry in the particle distributions in 
the hard-transverse jets around their axes. It is concerned with 
the large difference between the transverse and longitudinal (with 
respect to the beam axis) momentum components of the hardly 
scattered partons. This increasing of the longitudinal momentum 
component takes place due to the interaction of the transverse jet 
quark with gluon fields in the longitudinally oriented flux tubes. 
The value of the additional longitudinal momentum of the jet quark 
after crossing the flux tubes (<AP> m 1.7 GeV) is sufficiently 
large to be observable in experiments. Because this asymmetry is 
produced by the axially symmetrical field fluctuations, it should be 
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absent in both the cases of hadronic matter and QGP and may 
serve as a backgroundless signature of the flux tube stage in heavy 
ion collisions. 
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